Predictions from a stochastic polymer model for the MinDE dynamics in E.coli 
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The spatiotemporal oscillations of the Min proteins in the bacterium Escherichia coli play an 
important role in cell division. A number of different models have been proposed to explain the 
dynamics from the underlying biochemistry. Here, we extend a previously described discrete polymer 
model from a deterministic to a stochastic formulation. We express the stochastic evolution of the 
oscillatory system as a map from the probability distribution of maximum polymer length in one 
period of the oscillation to the probability distribution of maximum polymer length half a period 
later and solve for the fixed point of the map with a combined analytical and numerical technique. 
This solution gives a theoretical prediction of the distributions of both lengths of the polar MinD 
zones and periods of oscillations - both of which are experimentally measurable. The model provides 
an interesting example of a stochastic hybrid system that is, in some limits, analytically tractable. 

PACS numbers: 87. 16. A- Theory, modeling, and simulations 87.17.Ee Growth and division 87.10.Mn 
Stochastic modeling 

Keywords: Min proteins, polymer dynamics, stochastic oscillations, hybrid dynamical systems 



I. INTRODUCTION 

The spatiotemporal oscillations of the Min proteins 
in the bacterium Escherichia coli have been well stud- 
ied both experimentally and theoretically (see \l] for a 
recent review). Whereas many of the proposed mod- 
els @, H, 0, H, H, 0> Hi describe the oscillations as a 
self-organised emergent pro p erty of a reaction-diffusion 
system, only a few 0, TToL [Til [HJ address the poly- 
mer nature of the relevant proteins that is evident in 
both in vitro [lH, 03] and in vivo [HI, [l6[ experiments. 
Cytrynbaum and Marshall [ll| described the oscillatory 
dynamics of the Min proteins solely in terms of sim- 
ple polymer assembly and disassembly dynamics coupled 
with concentration-dependent conditions for switching 
between these two states. Much of the experimentally 
observed behaviour of wildtype cells can be explained 
within this framework as can a number of mutant-study 
observations. Here we analyze a stochastic version of 
the model introduced in [11] and build a set of quantita- 
tive predictions that can be used to evaluate the model 
against data and results of other models. 

Our stochastic model consists of a set of four interact- 
ing linear polymers, a pair of MinD and MinE polymers 
at either pole of the cell, each of which can be in either a 
growing or a shrinking state. For any fixed combination 
of states for the four polymers, the dynamics are deter- 
ministic and described by a system of ordinary differen- 
tial equations for the lengths of the polymers. The full 
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state of the system is thus determined by four discrete 
state variables (growing/shrinking) and four continuous 
variables (polymer lengths). Stochastic transitions be- 
tween the discrete states are dependent on the cytoplas- 
mic concentrations of the Min proteins and so indirectly 
on the polymer lengths. 

Cytrynbaum and Marshall [11] analysed a determinis- 
tic limit of infinitely high cooperativity, the solutions of 
which can be expressed easily in the form of fixed points 
of a one-dimensional map which we describe in the first 
subsection of Sec. [IVJ Formulating the model in this de- 
terministic limit as a map is useful in understanding our 
approach to the full stochastic model which we focus on 
throughout the rest of Sec. IIVI As with the determin- 
istic model, the stochastic model can be reduced to a 
map, in this case one that takes the probability distri- 
bution for the maximum polymer length during a given 
period of the oscillation to the probability distribution 
for the same quantity half a period later. Our results 
include analysis of this probability distribution map and 
the calculation of its fixed point. In addition, from the 
analytical theory and numerical simulations, we calculate 
the distributions of other properties that are easily mea- 
sured experimentally. This enables us to narrow down 
the parameter regime in which the results of our model 
agree with experimental observations. In particular, a 
high cooperativity [n ~ 6) in the nucleation of the MinE 
polymer is required. The cooperativity in nucleation of 
the MinD polymer turns out to be less crucial (n w 3 is 
sufficient). Also, over a wide range of parameter values, 
two stable solutions exist, one in which the Min proteins 
are entirely in the cytosol and one in which polymers 
form at either pole in an oscillatory manner. Stochastic 
transitions between these bistable states are more or less 
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likely depending on parameter values. In the discussion, 
we put this finding' into the context of recent experimen- 
tal observations [17|, and compare our results to other 
modelling studies in the literature. 



II. BIOLOGY 

It is well established experimentally that two main pro- 
cesses control the position at which the rod-shaped bac- 
terium E.coli divides. Nucleoid occlusion [1 81 ] prevents 
division at sites within close proximity of the nucleoid, 
which - at the relevant time after DNA replication - is 
everywhere in the cell except at the middle and near the 
two poles. The latter two sites are ruled out by the coop- 
erating action of a group of three proteins, the Min pro- 
teins MinC, MinD, and MinE. In E.coli, a spatiotemporal 
oscillatory pattern formed by these proteins restricts di- 
vision to the middle of the cell. 

MinD is an ATPase that, in its ATP-bound form, binds 
to the inner cell membrane [19j. MinE is found to acti- 
vate the ATPase activity of membrane-bound MinD and 
thereby removes MinD from the membrane [2(j. MinC 
co-localises to membrane-bound MinD [2l| and is known 
to prevent assembly of FtsZ, one of the important players 
in forming the apparatus that constricts the cell during 
division [22[ • Fluorescent labelling of MinD showed a pre- 
ferred localisation to the polar regions of the membrane 
and away from midcell [231 ] . This localisation appears to 
be the result of spatiotemporal oscillatory dynamics pro- 
duced by the membrane-dependent interaction of MinD 
and MinE (for a review see [l|]). The still- uncertain mech- 
anism by which this oscillation occurs is the motivation 
for our work. 

The interaction of MinD and MinE in the presence of 
ATP and lipid membranes has been studied in vitro and 
MinD was found to accumulate into polymers and fibre 
bundles above certain concentrations [14| . In vivo, the 
fluorescently labelled Min proteins appear to be organ- 
ised into helical structures on the inner wall of the cell 
membrane (TBI . [iH ] . The qualitative agreement between 
the in vitro and in vivo observation suggests that the Min 
proteins organise in the form of polymers in the cell, an 
assumption on which we base our theoretical model. 



III. MODEL - A HYBRID DYNAMICAL 
SYSTEM 

We assume that both MinD and MinE aggregate only 
in the form of polymers on the inner side of the cell's 
membrane. The Min oscillation results from an interplay 
between two MinD- and two MinE-polymers, one pair on 
each side of the cell. 

MinD monomers [13] from the cytosol can start a poly- 
mer (i.e. nucleate) at one of the nucleation sites that are 
assumed to be positioned at each pole of the cell, an idea 
supported by recent experiments [24l l25l . |26| . The proba- 



bility of nucleation on an empty site is proportional to the 
cytosolic MinD concentration raised to the power n lmc . 
We assume that each nucleation site, when occupied by 
a polymer, is incapable of nucleating a second polymer. 
The MinD polymer then elongates towards midcell at a 
rate proportional to the cytosolic MinD concentration. 
A MinE polymer can nucleate at the growing tip of the 
MinD polymer with a probability proportional to the cy- 
tosolic MinE concentration raised to the power n cap . It 
then grows backwards on top of the MinD polymer with 
a rate proportional to the MinE concentration. By in- 
ducing hydrolysis, the MinE subunits destabilise the un- 
derlying MinD subunits which disassemble from the tip, 
releasing both types of subunits into the cytosol. 

The geometry of our model cell is that of a cylin- 
der with fixed length L and diameter 2r. Reported 
diffusion coefficients for the Min proteins are around 
10 /Ltm 2 s _1 (2?| so in a cell of length 2 — 3 /im and with 
characteristic reaction times on the order of seconds, the 
cytosolic concentrations of MinD and MinE are essen- 
tially uniform throughout the cell. Also, the time scale 
of ADP-ATP exchange in cytosolic MinD is assumed to 
be fast compared to the oscillatory dynamics [Til ] . 

The equations governing the dynamics of the polymers 
are distinct for the different discrete states between which 
the system jumps. The system describing the behaviour 
therefore is a hybrid dynamical system - a combination 
of continuous and discrete dynamics (28[. In our case, the 
continuous variables are four polymer-lengths if, and if, 
whose dynamics are determined by the values of the re- 
spective discrete state variables Sfj r and Sfj r . In the fol- 
lowing, we will provide the model equations and rephrase 
the model introduced in [Tl[ using a slightly different no- 
tation. 



A. Polymer dynamics 

For each of the four polymers (two MinD, two MinE), 
we track the projection of the polymers onto the long axis 
(x) of the cell (Fig. [J). The four variables if, if , if and 
if describe these projected lengths for the polymers at- 
tached to the left and right pole of the cell, respectively. 
The relation between the full arc length of a helical poly- 
mer and its projection is given by 1/ cos 9, where 9 is the 
pitch of the helix. The parameter 7 = d cos 9 is used to 
convert between the projected length I of a polymer and 
the number of monomers it is made of (with monomer 
size d). 

1. MinD 

MinD polymerisation can start at either of two nucle- 
ation sites located at the two poles of the cell, i.e. at the 
positions x = and x = L. We assume that one end 
of the MinD-polymers is fixed to one of these stationary 
nucleation sites, whereas at the other end (the 'tip'), the 
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polymer can elongate or shorten. Growth and shrinkage 
are governed by the equations: 



(1) 



fc D (S' ; I y r ) represents either constant disassembly or first 
order assembly and depends on the discrete state variable 
Sy r of the polymer in question: 



on 



if< = 







if s. 



l/r 



(D-polymer growing) 
(D-polymer shrinking) 
1 (no D-polymer) 

(2) 

The dynamics of switching for the discrete state vari- 
ables S® r , between the three states —1,0 and 1, will be 
explained in the next subsection. 

The cytosolic concentration cd of MinD is determined 
by conservation of monomers: 



(I? + 1?) 



CD — CDto 



(3) 



where V = Trr 2 L is the volume of the cell (in /im 3 ), A ~ 
/^M/im 3 converts between particles per /Ltm 3 and /j,M, 
and cd, to is the total concentration of MinD monomers. 

The largest possible extension of a D-polymer is 
reached when all MinD is bound in one polymer. This is 
the case at 



T 



-CD, to- 



(4) 



For high total MinD concentrations, Z max can come close 
to L, i.e. the D-polymer would cover the whole cell from 
pole to pole. To avoid further assumptions on what hap- 
pens if a polymer hits the opposite cell wall, we restrict 
ourselves here to total MinD concentrations that make 
these events very unlikely or impossible. With the pa- 
rameters from Tab. [U this means we consider maximal 
total MinD concentrations of around 5 /jM. In the simu- 
lations, the polymer simply stops growing in the unlikely 
case that it reaches the opposite cell wall. 



2. MinE 

We assume that the MinE polymer nucleates on the 
tip of the MinD polymer and grows on top of it towards 
the pole. The differential equation describing the pro- 
jected length of the MinE polymer has the same simple 
structure as the one for the MinD polymer: 



dt 



-I. 



dt 



(5) 



The same conversion factor 7 is used since we assume 
that MinE monomers bind to MinD monomers of the 



helix one-to-one. The MinE-polymer always starts grow- 
ing from the non-polar tip of the MinD-polymer (Fig.[T]). 
The non-polar end of the MinDE-polymer falls off the 
membrane (with a slower speed than E-elongation) and 
disassembles. The state variable S , , E / „ defines the value of 
the growth rate as follows: 



l/r 



fc on C E 



k E (Sf /r ) 



Voir 



if S* E 

11 S l/r 

if ^ E 
11 °l/r 



1 (E-polymer growing) 
(E-polymer reached 
cell wall (lf /r = lf /r )) 



if Sf /r = -1 (no E-polymer) 

(6) 

As for MinD, we assume the total number of MinE 
monomers to be constant: 



_A_ 

7V 



(If + If) + c E = c E , to . (7) 



B. Switching 

We consider a stochastic description of switching be- 
tween the three possible discrete states of the state 
variables Sf/ r E where the probability of switching de- 
pends on the cytosolic concentrations of MinD/E. To 
capture the cooperative nature of the initiation of a 
polymer (e.g. [29|, HoJ), we assume that the instanta- 
neous rates of nucleation and capping are proportional 
to a power of the respective cytosolic concentrations: 
-M7r(*) = fc nucC D nuc (£) and A^(f) = fc cap c E cap (t). n nuc 
and n cap are the cooperativities for the nucleation and 
the capping events, respectively. A D-polymer starts 
growing out of an empty nucleation site with probability 
Pnuc(i) = A nU c(i)d£ during the time interval t..(t + dt) 
and a growing D-polymer gets capped (i.e. an E-polymer 
nucleates at its tip) with probability p C a, P {t) = A cap (i) dt 
during the time interval t..(t + dt). Our main results 
involve calculating analytical expressions for the proba- 
bility distributions when these switches happen. 

As a simplified model, we first consider the limit in 
which the cooperativities go to infinity. More specifi- 
cally, we define the nucleation probability as A™ c (t) = 

^nuc( c D(£)/cD,th)" nuc - In the limit n nuc -> 00, the 
stochastic switching becomes deterministic with the nu- 
cleation event occurring as soon as cd reaches the nucle- 
ation threshold co,th- Capping is treated similarly. 

Note that this deterministic limit of our model is 
fundamentally different from the deterministic model of 
Drew et al. [9y in that Drew et al. consider the mean 
field behavior of a population of filaments that switch 
between growing and shrinking states at average rates 
whereas we have individual filaments that switch at spe- 
cific concentration-determined times. The difference be- 
tween these two models is that the mean field model 
fails to admit oscillations without further biochemical 
assumptions (length-dependent growth speed), whereas 
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the individual-filament model has an oscillatory solution 
over a large range of parameter values. 



C. Model summary 

Eqs. H])-® together with the above-mentioned 
switching dynamics of the discrete state variables Sf/® 
represent a hybrid dynamical system [28]. In the rest of 
this article, we analyse this system and present both ana- 
lytical as well as numerical results that can be interpreted 
with regard to experimentally obtainable data. 

A single MinD-polymer 'life span' includes nucleation, 
growth, capping (nucleation of the MinE-polymer), and 
disassembly as described for the deterministic case in 
Tab.U 

The two poles can undergo alternating events of this 
type, which then constitutes an oscillatory solution. Such 
an oscillatory solution progresses in the following man- 
ner (cf. Fig. Q] and Fig. [3 for notation). A MinD polymer 
capped (at time tf) by a MinE polymer disassembles at 
one pole while the nucleation site at the other pole is 
either occupied by a disassembling MinDE polymer or 
remains empty. As the MinE polymer disassembles at 
the tip, it maintains a steady length by growing at its 
other end (treadmilling) . Throughout this process, cy- 
tosolic MinD concentration increases thereby increasing 
the probability of MinD-polymer nucleation at the empty 
pole but cytosolic MinE concentration remains low. Once 
nucleation occurs (at time t%), the nascent polymer grows 
towards midcell, depleting the cytosolic MinD pool. Be- 
cause of the ordering of critical concentrations for nucle- 
ation and elongation, the existing MinE polymer elon- 
gates in preference to nucleation of a new MinE polymer 
capping the nascent MinD polymer. When the original 
MinD polymer reaches the same length as the steady- 
state-treadmilling MinE length, treadmilling is no longer 
possible and the cytosolic MinE concentration begins to 
rise. This raises the probability of capping the nascent 
MinD polymer (at time f§). After capping, the system 
is back to the state we began describing but with the 
poles reversed. The first MinDE polymer completely dis- 
appears at time if. 

For a specific subclass of the described oscillatory so- 
lution we are able to derive an analytical description for 
relevant probability distributions. This subclass we call 
regular oscillations and it is defined by t\ < <t\< tf. 
A section of it is shown in Fig. [5J This definition essen- 
tially means that the growing phase (Sfj r = 1) of one 
polymer falls completely within the shrinking phase of 
the other (S^ fl = 0). 

Table [II] gives an overview and numerical values for 
the parameters we used in this model. Most of them are 
taken in ranges reported in the experimental literature. 
Some are adjusted to values that lead to reasonable re- 
sults of the model. 



IV. RESULTS 

A hybrid dynamical system must be solved piecewise 
and care must be taken at points of discontinuity which, 
in this case, occur each time a polymer switches state. 
Within the appropriate parameter ranges and for a reg- 
ular oscillation (as defined above), we need only consider 
the progression through three of the six possible com- 
binations of discrete states of (Sf,S®) which occur in 
the sequence (0,-1) — > (0,1) — > (0,0) during one half- 
period. The evolution of the continuous variable, for any 
given discrete state, requires a solution for disassembly 
and assembly of polymers. The solution of Eq. [1] for a 
shrinking D-polymer [Sy r — 0) is 



l?/r(t) = l?/r(t!/r)-lk oS (t-tf /r ) ) 



(8) 



where tf, r is the time of the most recent capping of the 
left or right D-polymer, respectively. 

The solution of Eq. [T] for a growing D-polymer is de- 
pendent on the dynamic cytosolic concentration. During 
regular oscillations, a growing D-polymer only appears 
while the other D-polymer is disassembling (i.e., discrete 
state (Sjj r , = (0, 1)). Assuming the polymer on the 
right is decaying, the cytosolic MinD concentration fol- 
lows c D (t) = c D , to - ±(l?(t) + l?(t c r ) - jk oS {t - O). 
Substituting this into Eq. [T] and solving the resulting 
equation with the initial condition £P(if) = (t™ is the 
time at which the left polymer nucleates), one obtains 
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7 V 



CD, to 



"■on 



7 y 



■ exp 



CD.to - Td" 

' on / , ,n \ 
■ — {f-ti ) 



l?(t C r )+jk oS (t-t C r ) 

(9) 



Similar solutions can be found for the dynamics of the 
E-polymers. However, for the analytical treatment in 
this article, we restrict ourselves to the limit of fast fi- 
ring formation: We assume that an E-polymer attains its 
steady state length right after its nucleation (i.e. after 
capping of the D-polymer). Equating the rate of polymer 
decay at the medial end of the MinDE polymer with the 
growth rate of the MinE polymer gives the steady state 

E-polymer length l E<ss = ^ (c E ,to - |ff )• With this 
simplification, the model reduces to a hybrid dynamical 
system with two continuous (ly r ) and two discrete (Sfj r ) 
variables. It is important to note that we make use of 
this approximation only for obtaining analytical results. 
In the numerical simulations we always model the full 
system with explicit E-polymer dynamics. 

For the rest of this article, we will drop the 'D' super- 
script and denote the length of the D-polymer by I. 



5 





deterministic switching condition 


5 D 


5 E 


(1) MinD polymer nucleates (nucleation event) 


Cd(£) > CD,th 


-1 -»• 1 




(2) MinE polymer nucleates (capping event) 


c B (t) > c E ,th 


1 -> 


-1 -> 1 


(3) MinE polymer reaches cell wall 


l B = Z D 




1 -> 


(4) MinDE polymer reaches cell wall 


P =0 


0-v-l 


-» -1 



TABLE I: The series of states a MinD polymer on one pole goes through during one 'life span'. 



' ~~ l Sf /r : / I 0/0 -1/0 1/0 1/0 

Sfy. / -1 0/1 -1/1 -1/1 -1/0 

FIG. 1: Scheme of the oscillations produced by the model. Shown are the two polymer helices (red - MinD, blue - MinE) on 
the left and the right side of the cell at different time points (time is increasing to the right). The values of the discrete state 
variables Sf,Sf, S? and Sf are given under the figures. 



par. 


value 


unit 


description 


d 


2.5 


nm 


increase in polymer length 
by addition of a single 
monomer (see footnote on 
page OH [3 





1.4 


rad 


pitch of the helical poly- 
mer [l5[ 


CD, to 


4 


^M 


total MinD monomer cone. 


CE.to 


1.5 


/iM 


total MinE monomer cone. 


feoff 


80 


s- 1 


depolymerisation rate of the 
MinDE polymer 


A- D 

"'on 


150 


(fM ■ s)- 1 


polymerisation rate of MinD 
on the membrane 


/,- E 

"'on 


320 


(fiM ■ s)- 1 


polymerisation rate of MinE 
on the MinD polymer 


CD,th 


2.5 


fiM 


threshold cone, for the nu- 
cleation of a MinD polymer 


CE.th 


1.25 


AtM 


threshold cone, for the nu- 
cleation of a MinE polymer 


L 


3 


fim 


length of cell 


f 


0.5 


fim 


radius of cell 


fenuc 


0.015 


s-VM" nnuc 


rate constant of nucleation 


fecap 


0.15 




rate constant of capping 


n nuc / 3-6 




cooperativity of nucle- 


^cap 






ation/capping 



TABLE II: List of parameters used throughout this article. 
The concentrations CD,to and CE.to are consistent with values 
used in other modelling papers and are supported by exper- 
iment [3j|. The rates governing polymer growth and decay 
are estimated from [32l ] . L and r are typical values seen in ex- 
periments and the other parameters are chosen such that the 
model produces reasonable results. A detailed discussion of 
the numerical values of some of the parameters can be found 
in [H. 



T ™ T* 




IT 

FIG. 2: A section of the regular oscillation pattern with 
stochastic switching to explain our notation. Shown is the 
temporal dynamics of the tip-positions (x) of the D-polymers 
as well as the value of the discrete state variables SY/ r (shaded 
areas). We use capital T for time differences and little t for 
time points. Even indices refer to the polymer anchored at 
the right pole of the cell and odd indices refer to the one 
anchored at the left pole. 



A. Solution of the deterministic model 

As a basis for further discussions we briefly present 
the solution to the simplest version of the model de- 
scribed in the preceding section and previously addressed 
by Cytrynbaum and Marshall [ll| . We consider the case 
of deterministic switching (see Subsec. IIII Bp and fast fi- 
ring formation. The length of the MinD polymer at cap- 
ping (the amplitude of the oscillation) on one side can 
be expressed as a function of the capping length at the 
preceding capping on the other side: = f(lf) (cf. 
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2. An unstable oscillation (the map intersects 
identity line with a slope larger than one). 



the 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

(a) varying c Djto 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 










^ cs.to ^''1.3 (iM 




/ 1.5 















0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

(b) varying c E ,to 

FIG. 3: The dynamics of the deterministic version of the sim- 
plest model (infinitely fast E-ring formation) depicted in form 
of a discrete map. It displays the length of the MinD poly- 
mer relative to the cell length L as a function of this lengths 
at the previous capping, i.e. the amplitude of the oscillation 
(see App.fX]for the equations). For intuitive analysis of the 
map, the identity line is added. Except for crj,to and CE,to, 
standard parameters from Tab. lU were used. 



Fig. [2]). Thus, the problem of finding a periodic solution 
to the hybrid system is reduced to finding a fixed point 
of a one-dimensional map. Only in this simple version 
of the model such a one-dimensional map can be found, 
because only then the state of the system is completely 
determined by the length of only one of the polymers. 

The calculation and the equations for the map are 
given in App. [AJ In Fig. [3] the map is plotted for varying 
total concentrations of MinD and MinE. 

The intersection of the map with = If shows three 
fixed points which in terms of the hybrid system corre- 
spond to: 

1. A stable cytosolic solution: There are no polymers 
(l c = 0) and all the Min proteins are in the cytosol 
as monomers. 



3. Stable (slope is equal to zero) oscillations with con- 
stant amplitude (given in Eq. IA.3|) . 

Depending on the initial conditions (k(t = 0)), the sys- 
tem converges to one of the two stable states, i.e. each 
parameter set that allows for an oscillatory solution also 
includes a solution where no polymer exists (the cytosolic 
solution) . 

As can be seen from Fig. [3J the region of initial condi- 
tions that lead to a stable oscillatory solution is biggest 
when cd, to is large and CE,to is close to CE,th- If CD.to is too 
small or CE,to is too high, the map does not intersect with 
the identity line anymore and only the cytosolic solution 
remains. This corresponds qualitatively to the condition 
for the existence of oscillations as derived in [ll| for a 
simplified version of the model. 

In the stochastic version of the model, low probability 
switching events can, under certain circumstances, lead 
to a transition between episodes of oscillations (with rel- 
atively constant amplitude) and almost purely cytosolic 
states. We will investigate this in the following sections 
and will refer to Fig. [3J as the limiting case of infinite 
cooperativity. 



B. Probability distributions for experimentally 
measurable quantities 

More realistic than the deterministic-threshold limit 
for nucleation and capping of the MinD polymer is 
stochastic switching that models the random nature of 
the underlying chemical reactions. From now on, we will 
use the stochastic rules for switching between the dif- 
ferent states of polymers as introduced in Subsec. IIII Bl 
In this section, we will characterise the changes such a 
stochastic switching rule introduces to the system. We 
will focus on two measures to characterise the robust- 
ness of the oscillatory solution of the stochastic model: 
the distribution of oscillation amplitudes and periods and 
the possibility of skipping beats. Both of these measures 
should be easy to obtain from experiments. 

As long as the system undergoes regular oscillations, 
conditional probability distributions can be derived for 
the times at which nucleation and capping occur. Sup- 
pose the polymer on the right disappears at time t d 
(Fig. and, at that moment, the polymer on the left 
is shrinking and has length l d . Given this, we denote the 
probability that a new polymer nucleates on the right at 
time t by P nVLC (t\l d ). Thus, P nuc (t\l d ) dt is the probabil- 
ity that nucleation on the right side happens in the time 
interval t..[t + dt) where t > t d . Similarly, if t n is the 
time of nucleation on the right and the polymer on the 
left is still shrinking and has length l n , then P cap {t\l n ) dt 
describes the probability for the capping of the growing 
polymer to happen in the time interval t..(t + dt) after 
nucleation (t > t n ). Analytical expressions for these two 
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probability distributions are computed in Apps.lBlandlCl 
and plotted in Fig. [TT] 

Using the two conditional probability distributions for 
nucleation and capping, analytical expressions for the 
steady state probability distribution of polymer lengths 
at capping and other relevant quantities can be derived. 
Similar to the description of the deterministic system in 
terms of a map (Subsec. IIV A [I . we derive an integral re- 
lationship that maps the probability distribution for the 
polymer length l n (Fig. [2]) onto a new probability distri- 
bution of the same length half a period later. The deriva- 
tion of this map P(lf) = F(P(lf_ 1 )) and the expression 
for the map itself (Eq. ID.6[) is presented in App. |Dj The 
expression for the operator F that relates two consec- 
utive probability distributions involves complicated in- 
tegrals and no general closed expression can be found. 
However, through numerical integration, an approxima- 
tion to the steady state probability distribution can be 
computed iteratively. Typically, three or four iterates of 
the map F (Eq. ID.6[) are adequate for convergence with 
reasonable accuracy. 

Only in the case of regular oscillations can the two 
stochastic processes, nucleation and capping, be sepa- 
rated and the probabilistic map (Eq. ID.6|) be derived. 
Crucial for obtaining regular oscillations are high coop- 
erativities in both capping and nucleation. 

In the following, we present probability distributions 
for the amplitude and period of the oscillations obtained 
from numerical simulations [35| . We compare these re- 
sults to analytical results that we obtain by iterating the 
probabilistic map. 



1. Distribution of oscillation amplitude and dependence on 
cooperativities 

Fig. 0] shows the results of long simulation runs of the 
full MinDE polymer model with stochastic switching. 
Plotted is the distribution of MinD-polymer lengths at 
capping (the amplitude of the oscillation) and the time 
between consecutive capping events on the same side (the 
period). 

In experimental studies (both qualitative and quanti- 
tative [3llp. large deviations in the oscillation amplitude 
and period are seldomly reported. For our model to re- 
produce this narrow distribution, high cooperativities in 
both nucleation and capping are crucial. Fig. [4] shows 
that our model is especially sensitive to the cooperativ- 
ity in capping: n nuc = 3 still leads to a fairly narrow 



distribution, whereas both curves for n c 



3 have big 



contributions at unusually short capping lengths. 



Fig.Oshows the same simulation results as of Fig. 4(a) 



together with the analytical result from the numerical 
iteration of the probabilistic map fEqs. IdTBI and |P77|) . 

The plots confirm the applicability of our analytical 
solution as discussed above. For regular oscillations, the 
analytical result agrees well with the results from sim- 
ulations. When cooperativities decrease, the increased 
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FIG. 4: Numerically obtained distributions for oscillation am- 
plitude (polymer lengths at capping) and periods (time be- 
tween consecutive cappings on the same side) with varying 
cooperativity of nucleation (MinD) and capping. Standard 
parameters as of Tab. [TT] were used. 



probability at small polymer lengths is not accounted for 
by the analytical result. These are the events that are 
too far away from the stable oscillation point in Fig. [3] 
and violate the regular oscillation assumption needed for 
the calculations in App. |Dj Since both the numerical 
and the analytical distributions are normalised, the lat- 
ter one shows higher values in the region corresponding 
to regular oscillations. The broader capping probability 
distribution in the case of lower n cap (cf. also Fig. 11(b) ) 
can lead to capping at short polymer lengths, which even- 
tually can lead to a transition to a purely cytosolic state 
(see Subsec. HVjJ. 
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FIG. 5: Comparison of the analytical result (line) with simu- 
lations (points - same data as in Fig. |4(a) |. For the analyti- 
cal result, the probabilistic map (Eq. ID.6"|l was applied three 
consecutive times, starting with a uniform distribution. This 
steady state distribution in l n was then transformed into a 
distribution in l c by applying Eq. ID.7I 



Distribution of oscillation periods and 
total concentrations 



Another quantity that can be easily obtained from ex- 
periment is the period of the oscillations. In Fig. [51 
we plot the probability distributions of the period T for 
different total concentrations of MinD and MinE as ob- 
tained from simulations. 



At higher concentrations of MinD (cd,i 



> 



3.75 (M) 



and intermediate concentrations of MinE (cE,to ~ 
1.5 /iM), the distributions arc sharply peaked around a 
single value of T (see Fig. 6(a) ). This parameter regime 



leads to regular oscillations. The amplitude of the oscil- 
lations decreases monotonically with c p^o, reflected in 
the leftward shift of the peak in Fig. 



6(a) 



consistent 

with the analogous feature in the deterministic version 
of the model (Fig. [3]). The closer the amplitude comes 
to the unstable fixed point in Fig. [3J the more likely a 
D-polymer will be capped at a small length. This early 
capping at one pole leads to an earlier rise in MinE so 
that the subsequent capping event of the other pole oc- 
curs earlier as well. This positive feedback can drive the 
cell into the predominantly cytosolic state which appears 
in the stochastic version of the model as a second peak 
close to T = (Fig. |6(a)j c D ,to = 3.25 (M). 

A similar distribution is obtained if the total MinE 
concentration is chosen too high (shown in Fig. |6(b)| 
CE,to = 1-6 /iM) with capping of the D-polymers occur- 
ring quite early. For low concentrations of MinE, the 
instantaneous capping rate is reduced and capping of the 
D-polymer can occur late (cf. also Fig. 11(b) I. In ex- 




100 

T[s] 
(a) varying c Djto 



200 




150 200 250 300 350 

T[s] 

(b) varying c E ,to 

FIG. 6: Distributions of oscillation periods T for different to- 
tal concentrations CD.to (upper panel) and CB,to (lower panel). 
The data was obtained from simulations where the time be- 
tween two consecutive cappings on the same side is considered 
to be T (ignoring the dynamics on the other side). Both coop- 
erativities (for nucleation and capping) are 6, CE.to = 1.5 /iM 
(upper panel), and co,to = 4/tM (lower panel). 



treme cases (cE.to = 0.8 /iM) the D-polymer can reach 
maximum length i max (Eq. 0]) . Accordingly, the distribu- 
tion of periods is shifted towards higher T and tails out 
slowly. 

Using the steady state probability distribution for l n , 
an analytical integral expression for the probability dis- 
tribution of the period T can be derived for the case of 
regular oscillations (see App.[E|. Fig. [7] shows a compar- 
ison of this analytical result with simulation data. Since 
the analytical description only contains the regular os- 
cillations, it fails to produce the peaks at short periods 
in Fig. 7(a) and |7(d)] Also, low total concentrations of 
MinE (e.g., CE.to ^5 1.2 /iM) lead to an irregular oscil- 
lation regime by violating the scheme shown in Fig. [2] 
in that capping of a new polymer is likely to occur after 
complete disassembly of the existing polymer at the other 



9 



0.04 






a, 0.02 




a. 


20 40 60 80 100 

T[b] 

(a) 3.25 aiM/1.5 fj,M 


0.04 












a. 002 




1 




a, 



0.04 



0.02 




25 50 75 100 125 

T[s] 

(b) 4 /jM/1.5 jtM 



0.04 



&, 0.02 




100 200 300 

T[b] 

(c) 4/xM/0.8/iM 



25 50 75 100 125 

T[b] 

(d) 4 ^M/1.6 fiM 



FIG. 7: Comparison of the analytical result (line) with sim- 
ulations (points). For the analytical result, the probabilistic 
map (Eq. ID.6[) was applied three consecutive times, starting 
with a uniform distribution. This steady state distribution in 
l n was then transformed into a distribution in the period T 
by applying Eqs. IE. II and IE. 21 The simulation data in (a), 
(c), and (d) is from Fig. [6] the data in (b) from Fig. [4] The 
numbers in the subcaptions are CD,to/cB,to- 



pole (through t\ > tf_ x ). Accordingly, the anal ytica l so- 
lution fails to give a good approximation (Fi g. |7(c)| . In 
Fig. |7(b)[ the standard parameter set (Tab. III[) is used 
and analytical and numerical results agree. Fig. |7(b)| 
and |5 (a)] represent the parameter values that provide the 
closest match to the experimental observations of sus- 
tained and regular oscillations. 

To illustrate the range of validity of our analytical cal- 
culation, we compare the means of the numerically and 
analytically obtained probability distributions in Fig. [8] 

The deviations in mean period in Fig. [S] reiterate that 
the analytic expression is valid provided MinD concen- 
tration is high and MinE concentration is intermediate, 
ensuring regular oscillations. For small CD,to, the analyt- 
ical result is always an overestimate of the numerical re- 
sult, since the analytical treatment does not capture the 
second peak in the distribution at short T (Fig, 
The same overestimate results for high CE,to (Fig 
whereas for low CE,to, the analytical calculation pr 
an underestimate of the period 



7(a) I. 
7[dJ|, 



rovides 
Here, the long tail 
of the numerically obtained distributions is not covered 



(Fig. 7(c)) 



Our results as shown in Fig. [8] agree qualitatively with 
experiments. Overexpression of MinD has been shown 
to increase the oscillation period 21( and overexpres- 
sion of MinE is known to disrupt the normal division 
placement 133], probably due to destroying the oscilla- 
tory pattern (which corresponds to our numerical solu- 
tion of T w s for high MinE concentration). In the 
deterministic version of the model, this corresponds to 
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FIG. 8: Mean oscillation period T for varying total concen- 
trations of the two Min proteins. Both the analytical re- 
sult (line) as well as results from simulations (points) are 
shown. Cooperativities for nucleation and capping are both 
6, CE,to = 1-5 /iM in the upper panel and CD,to = 4 /iM in the 
lower. 



the deterministic map not intersecting with the identity 
line anymore (Fig. 3(b)). 



3. Beat-skipping 

Stochastic switching (nucleation and capping) can also 
lead to occasional skipping of beats, another property 
that is easily quantified experimentally. Regular oscilla- 
tions are characterized by the asynchronous and alter- 
nating growth and disassembly of a polymer at each of 
the two poles of the cell. A skipped beat is defined as 
a deviation from regular oscillations in which alternation 
fails and a polymer reappears on the same pole before 
the appearance of a polymer on the opposite pole. As a 
criterion for our numerical analysis, we count the number 
of cases where two consecutive cappings of MinD poly- 
mers occur on the same side of the cell. The fraction of 
these cases with respect to all cappings is shown in Fig. |9] 
The figure shows numerically-obtained data for different 
cooperativities of nucleation and capping. 

As a skipped beat is one way of deviating from regular 
oscillations, the results in Fig.[9]recapitulate the observa- 
tions described earlier in the context of Figs. HI El and [8] 
The high fraction of irregular capping events for low to- 
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C. Bistability and stochastic transitions 
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FIG. 9: Numerical result for the regularity of the oscillations. 
Plotted is the fraction of skipped beats (see text) depending 
on both total concentrations of the Min proteins and coop- 
erativities used for nucleation and capping, respectively. In 
the upper panel, CE.to = 1-5 (jM, in the lower co.to = 4/xM. 
Simulations were run such that the total number of cappings 
are in the range of 10 s for the low values and 10 3 — 10 5 for 
the high values. The errorbars are smaller than the pointsize 
and the lines are guides to the eye. 



tal MinD concentrations as well as for high total MinE 
concentration correspond to the second peak in the prob- 
ability distributions of periods in Fig. [5] at short T. If the 
capping cooperativity is relatively low, there will always 
be a significant fraction of skipped beats, which corre- 
sponds to the non-zero tails of the distributions at short 
lengths for n cap = 3 in Fig. 3J For our standard param- 
eter set n nuc = 6,n cap = 6,c D , t o = 4^M, c E ,to = 1.5 /jM 
and for a range of concentrations around it, the fraction 
of skipped beats is well below the 1% figure. The strong 
deviation from our analytical result when CE,to is very 
small (Fig. 7(c) ) is a violation of the regular oscillation 
pattern that does not lead to skipped beats, i.e. even in 
this extreme case the alternating appearance of a poly- 
mer on the left and the right side is still conserved. 



As derived in Subsec. IIV Al the deterministic system 
is bistable for a large range of parameters. Under cer- 
tain conditions, stochastic nucleation and capping of 
polymers can lead to transitions between the two sta- 
ble states: oscillatory episodes lasting for tens of peri- 
ods and the cytosolic state. Fig. [10] shows examples of 
this behaviour for three different parameter sets. A clear 
distinction between regular oscillations and the cytosolic 
state is only possible for high cooperativity. This be- 
haviour can be explained qualitatively using the deter- 
ministic map described in Subsec. IIV Al 

During regular oscillations, a D-polymer is most likely 
capped at a length close to the peak of the respective con- 
ditional probability distribution (Fig. 11(b) in App. [Cjl . 
Due to the stochasticity of the capping process, there is a 
small probability of the D-polymer being capped before 
it reaches normal extension |36J. For high cooper ativi- 
ties, the deterministic map in Fig. [3] retains some validity 
and one can think of the capping events in the stochas- 
tic model as being a blurred version of the deterministic 
map. Rare early and late capping events correspond to 
events at the edge of the blurred region about each fixed 
point and can result in a transition from relatively stable 
oscillations to a state that is mostly cytosolic and back 
again. Since the events allowing for such transitions lie 
outside the regular oscillation pattern (Subsec. IlII C[) . the 
analytical approach in App. [D] does not provide the right 
tools to describe this effect. We restrict ourselves to a 
qualitative discussion of a limiting case. 

In order to get a better understanding of the stochastic 
transitions we consider the special case of deterministic 
nucleation of D-polymers and stochastic capping (37j. For 
high values of n cap , the map in Fig. [3] can be used as a 
rough guide for the stochastic transitions. In App. [0 we 
adjust this map to the specific case considered here. 

In Fig. [I0(d)1 we show three maps obtained as in 
App. |F) One representing our standard parameter set 
CE.to = 1.5 fM, fccap = 0.15 8~>M - 6 (cf . Fi g.|7fb)T ) and 
the two parameter sets used in Fig. 10(a) and 10(b) (with 
CD,th = 2.5 /iM). One can think of the stochastic system 
as producing a cloud of points in this map around the 
stable steady states. A transition occurs if an extreme 
capping event takes the system from a neighbourhood of 
the pre-transition fixed point over the unstable steady 
state into a neighbourhood of the other fixed point. For 
our standard parameter set, due to the large distance 
between stable and unstable fixed points, this is unlikely 
to happen. The system will stay in the vicinity of one 
of the fixed points with probability close to one. When 
the total E-concentration is increased (Fig. 10(a)), the 



map is shifted to the right and the distance an extreme 
capping event has to cover in order for a transition to 
occur, decreases. The same is valid for an increase in 
the capping parameter fc cap as observed for Fig. |10(b)| 
The map for this case is slightly shifted to the left com- 
pared to the former case which leads to a reduced transi- 
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FIG. 10: Episodes of oscillatory dynamics alternating with the cytosolic state. Plotted is the position (along the cell long axis 
x) of the two D-polymer tips over time. Note that, typically, the duration of one period in the oscillatory state is around 
T ~ 100 s. The map shown in (d) is used for a qualitative explanation in the text. Note, that contrary to our discussion in the 
text, the time series presented in (a)-(c) were obtained from the full stochastic model (nucleation and capping stochastic). 



tion rate. The main parameter controlling the extension 
of the cloud of points is the capping cooperativity n cap . 
When it decreases, the cloud becomes larger, rendering 
the map interpretation less meaningful (cf. Fig. rj| . More 
transitions occur and the states are not as well defined 



anymore (Fig. 10(c) I 



V. DISCUSSION 

The mechanisms underlying the Min oscillations in 
E.coli are still subject to much debate among modellers 
since none of the models proposed so far capture all the 
properties and all the phenotypes displayed in wildtype 
and mutant cells (see [ll| for a discussion). The most ob- 
vious qualitative features of the observed dynamics are 
explained in many models which means that being able 
to distinguish between these models will rely on careful 
quantitative characterization of both the models and the 
experimental data. The results presented here are a step 
in this direction. 

We generalized a recently published model to allow 
for more realistic comparison with data by introducing 
and analyzing stochasticity. From a mathematical view- 
point, this model provides an example of a stochastic 
hybrid dynamic system, whose solution can be found an- 
alytically in the case of regular oscillations. Probability 
distributions for easily measured quantities are provided, 
making the model testable. Comparing our results to ex- 
perimental data already available in the literature, the 
main findings of this paper are as follows. 

High cooperativity in the nucleation of the Min- 
polymers (at least a power of 3 for the MinD and 6 for 
MinE) are important. Lower cooperativities allow too 
many events that are inconsistent with observations. Re- 
ported Hill coefficients for the nucleation of MinD on a 
lipid membrane are around 2 [29| . 

The model is relatively stable to changes in MinD con- 



centration. MinD concentrations greater than 3.5 /iM 
ensure regular oscillations. MinE concentration, on the 
other hand, is more finely constrained and must fall be- 
tween 1.2 and 1.6 /jM. 

Our results offer an explanation for the variability in 
the presence of oscillations observed in experiments. A 
significant fraction of cells usually do not show oscilla- 
tions in experiments or go from an oscillatory state into a 
cytosolic state during the course of observation (e.g. [HI). 
Our model suggests an explanation for this in the form 
of bistability of an oscillating state and a purely cytosolic 
state. The model of Fange and Elf Q also reports bista- 
bility, whereas the models of Howard et al. Q and Meacci 
et al. demonstrate qualitatively distinct solutions (cy- 
tosolic or oscillating) as a function of parameters but not 
bistability. Transitions between the two states will be 
a challenge to observe experimentally. Evidence for a 
transition in each direction would be necessary, requir- 
ing long observation times. Photobleaching in the fluo- 
rescence microscopy studies limits the latter and other 
factors might influence the existence and quality of the 
observed oscillations (see, e.g. HH). 

In a recent experimental study, Downing et al. [TtJ ac- 
tively controlled the MinD dynamics in E.coli by chang- 
ing cationic concentrations in the surrounding medium. 
They were able to stop, distort and restart the oscilla- 
tions. Based on the bistability and stochastic transitions 
described here, we speculate that the ion concentration 
has an influence on the assembly dynamics of the MinE- 
polymer. For example, if increased ion concentration in- 
creases fc cap , the system might be driven into a regime 
where a stochastic transition into the cytosolic state be- 
comes more likely. On the other hand, decreasing fc cap 
can lead to a freeze or or halt of the oscillations with 
most of the MinD localised on one side of the cell. 

Most other models [!, [T3, EH produce a similar depen- 
dence of the oscillation period on the total Min concen- 
trations and a few [|| 0] report on the variability of the 
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period or, as in our case, the full probability distribu- 
tion P3. 

The biggest difference between the reaction-diffusion- 
type models and the polymer model described here is that 
in the polymer model the spatial pattern of the Min os- 
cillations is prescribed in the form of possible nucleation 
sites of MinD in the membrane. This assumption be- 
comes most important when trying to model the striped 
Min phenotypes observed in filamentous cells [23| . There 
is an ongoing debate in the experimental literature about 
the existence and possible nature of sites to which MinD 
binds preferentially. Recent work [24j has shown that 
specific anionic phospholipids colocalise with MinD at 
the poles and septum. An important as yet unanswered 
question is which of these two determines the localisation 
of the other or if there is some other factor upstream of 
both. 



APPENDIX A: THE DETERMINISTIC VERSION 
AS A ONE-DIMENSIONAL MAP 

In its simplest version (deterministic switching and fast 
E-ring formation, cf. beginning of Sec. IIV|) . our model 
allows for a fully analytical solution in terms of a discrete 
map. We choose to display the maximal length of the 
MinD polymer as a function of the temporally previous 
maximal length on the other side of the cell: l\ , x = f{lf) 
(cf. Pig. p. 

To derive the map, we consider the following case: At 
time t = 0, the DE-polymer on the left side (I) is decaying 
(Sf(0) = 0) and the D-polymer on the right side (r) is 
growing. Without loss of generality, we choose t — such 
that the DE-polymer on the left just reached the length 
h(0) — ^x"(cE,to — CE,th), which causes the cytosolic E- 
concentration to grow above CE,th- This causes the D- 
polymer on the right side to be capped: t c r — and 
therefore its state to be switched: S^! : 1 — > 0. We 
assume that the D-polymer on the right had a length 
Z r (0) = J£ at the time point of capping. For a short 
period of time, now both polymers are decaying: lu r (t) = 
h/r(Q)—jk st. At time tf = fc * - h(0), the left nucleation 
site becomes free for a new nucleation. 

We now have to distinguish three different cases: 

1. The newly nucleated polymer on the left side gets 
capped immediately if CE(if) > CE,th- This is the 
case if 



0; tf 



tf), as long as 
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< 2 — (cE,to - C E ,th)- 



(A.l) 



Under this condition, the cytosolic MinE concen- 
tration will never drop below CE,th- Any nucleating 
MinD polymer will therefore be capped immedi- 
ately, which evolves into the cytosolic state (i.e. no 
polymer: If — 0). 

2. For slightly larger i£, nucleation on the left side 
will happen right after the first polymer (left) com- 



pletely disassembled (k 

c D (if) > c D , th - Using t\ = ^{l c r ~ ^(ce.to - 
c E,th)] m Eq. EJleads to an expression for the length 
of the D-polymer at the next capping: 



CD, to 



CD. to 



fcog 
jUD 

fcoff 

jUD 



CE.to + CE.th 



" CE,to — CE.th 



• exp 



lk off V 



lr — 2-^-(CE,to 



c E,th, 



(A.2) 



3. If l c r is even bigger (1% > ^(e D , t h + c E ,to - c E ,th)), 
such that cd (if) < CD.th, the new polymer on the 
left will nucleate at tf = ^^[Z£-^(cD,to-CD,th)]- 
Its capping will then happen at the stable ampli- 
tude: 



CD, to 



CD.th 



l,D 

"■or 

fcoff 
jLD 

"<on 



CE.to + CE,th 



• exp 



ft on 
fcoff 



(CD.to " CD,th — CE.to + CE,th) 



(A.3) 



This is the amplitude of the stable oscillation and 
all consecutive cappings will happen at this ampli- 
tude, too. 

The piecewise map derived above is displayed for var- 
ious MinD and MinE concentrations in Fig. [31 



APPENDIX B: CONDITIONAL PROBABILITY 
DISTRIBUTION FUNCTION FOR NUCLEATION 
TIMES (REGULAR OSCILLATION) 

In the stochastic version of our model, the probability 
distribution function of the times one side of the cell is 
free of polymer can be derived analytically if the system 
is oscillating regularly (see Subsec. Ull C[) . This distribu- 
tion function is conditional on the length l d of the poly- 
mer present on the opposite side of the cell at the time 
point where the nucleation site becomes polymer free (t d 
in Fig. [2]). As described in Subsec. IIII Bl we assume the 
instantaneous probability for nucleation on a free nucle- 
ation site at time t to be A nuc (t) = fc n ucCQ nuc (i). 

The probability for nucleation between times t and 
t + dt is the probability of no nucleation until t and 
then nucleating in t..(t + dt): P nuc (t) dt = (1 - P(t n < 
t))X nuc (t) dt. P(t n < t) is the cumulative distribution 
function. Without loss of generality, we assume tf = 
(i.e., the polymer on the left side just decayed completely 
at time t = 0) and obtain 



-Pnuc(i) = A nuc (i)exp 



Kuc(t')dt' 



(B.l) 
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The cytosolic D-concentration follows Co(t) = Cn.to — 
^(lf ~ iK&t) for t < tf (on the other side). For t > tf, 
there would only be cytosolic MinD (co(t) = co,to)- We 
do not consider this case, since it is outside the regular 
oscillation pattern. 

Putting cd (t) from above into Eq. IB. 11 one obtains 



Pnuc(T f \l d ) 



■ exp 




Fig. [TT1 shows typical shapes of this probability distri- 
bution function (for a given l d ) for high cooperativity in 
nucleation and different concentrations of MinD. From 
Eq. IB. 21 it is obvious that the dependence of the prob- 
ability distribution on l d is (up to a scaling factor) the 
same as on CD.to- 



APPENDIX C: CONDITIONAL PROBABILITY 
DISTRIBUTION FUNCTION FOR CAPPING 
TIMES (REGULAR OSCILLATION) 

Equivalently to the preceding appendix, one can de- 
rive an expression for the conditional probability distri- 
bution function for the time of growth of a polymer. 
In Subscc. IIII Bl we introduced the instantaneous prob- 
ability for capping of a growing polymer at a time t: 
A(i) = k cap o^ aT, (t). 

We only consider the approximative case of fast E-ring 
formation (see Sec. II VI) . The MinDE-polymer on the 
right side is assumed to be decaying and has length Z™ 
at time t — (cf. Fig. [2]). Without loss of generality 
we assume that the D-polymer on the left side nucleates 

and starts growing at time tf = 0. We now have to PnQf) — k nuc ( cn.to 
distinguish three cases during the decay of the polymer 
on the right side (t d ' E is the time when the growing tip 
of the E-polymer reaches the cell wall): 

1. First, the concentration of MinE monomers in the 
cytosol is constant (if = ?e, ss => ce — ce,o = 

feoff 



2. For tf' E < t < tf, ce grows linearly: ce(£) = ce,o 
$k oS (t-tf> E ). 

3. After tf, ce is constant again: ce = cs.to- 



with Ie.ss = -^-(cE.to — fei? 1 )- The condition on l n only 
enters through the dependencies in t d ' E and t d . 

Fig. llll shows typical probability distribution functions 
for capping (for a given l n ) for high capping cooperativ- 
ity and different MinE concentrations. The dependence 
of the probability distribution on l n is (up to a scaling 
factor) the same as on CE,to for the important first and 
second part of the non-smooth Eq. 10. It 



APPENDIX D: A MAP FOR THE PROBABILITY 
DISTRIBUTION 

Combining the two conditional probability distribu- 
tions computed in Apps. iBlandlCl one can derive a map 
for the probability distribution of the amplitude during 
regular oscillations. Here, we derive this map in the fol- 
lowing form: P n (l\ l ) = F(P n (Z?_ x )), i.e. the probabil- 
ity distribution of lengths at nucleation (on the opposite 
side) can be derived as a function of the distribution at 
the previous nucleation (on this side - for notation see 
Fig. 02). 

We start with computing a relation between the two 
probability distributions P n {V^) an d PdQf )'■ 



' P (i?\it)p d (it) ai° 



(D.l) 



The conditional probability distribution p{lf\l d ) can be 
obtained from Eq. |R2]by replacing T f = -^Qt - l ?) 
and multiplying with (to ensure correct normalisa- 
tion under the variable transformation). For the limits 
of the integration one has to consider that If will always 
be smaller than if (see Fig. [2]). The upper limit is the 
maximally possible length of a D-polymer (i.e. all MinD 
is bound in one polymer (Eq. [4])). 
With that, we obtain: 



A 



7 y 1 



CD,to 



exp 



CD, to 



exp 



1 



Vk u 



+ 1 Afcoff 



Vk n 



1 Afc ff 



Pd(lf)dl d . 

(D.2) 



Putting these into the equivalent of Eq. IB. II gives the 
conditional probability distribution function for the time, 
at which the growing D-polymer on the left side gets 
capped and switches states (see Eq. IC.ll - we generalise 
to T c and drop the side-dependence; ce,o = |%F; A — 



In a similar fashion, Pd(lf) can be expressed as a func- 
tion of P„(J?_i): 



P d (lf) 



(D.3) 



Afc ofE ). Here, t d = ^l n and t d > E 



7feoff 



(l U ~ Ie,ss) 



The conditional probability distribution in Eq. ID. 31 can 
be obtained from Eq. IC.ll when replacing T c with the 
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p c , P (T c \n = < 



E,0 1 I 



fccap(CE,0 + A(T C -i d ' E )) 



,o + A{T - r^))'"^ 

(-fecapC^ B - [(CE,0 + A(T< - t^))- P + l _ 

"cap / , "cap ,d,E _ fc cap 

,ra P C E,to ex P ^ ^capC-jj g t A(n cap + 1) 

•[(c Ej o+A(t d -^ E ))—+ 1 - C ^ p+i ] - 



< T c < t" 



t d,E < T c <t d ^ ^ 




2 4 6 8 10 12 14 16 

T f [si 

(a) Pnuc{Tf\l d ), l d = 1.5 Mm 




T c [s] 

(b) Pca P (T c |r), r = 1.5 Mm 

FIG. 11: Typical probability distributions for the time free of 
polymer T (top) and the duration of growth T c (bottom), 
respectively. A number of curves are plotted for different total 
concentrations of MinD or MinE, respectively. The analyti- 
cal expressions are given in Eq. IB. 21 for the upper panel and 
Eg. IC. II for the lower one. Standard parameters as of Tab. ITT1 
are used and the probability distributions are given for the 
conditions l d = 1.5 /im and Z n = 1.5 fim, respectively. The 
cooperativity for nucleation or capping is 6. For comparison, 
the dotted lines show the probability distribution functions 



for a reduced cooperativity n n 



c / cap 



3. Note the broader 



distributions for lower cooperativity as well as the noticeable 
tail towards shorter T c in the lower panel for the dotted line. 



solution of the deterministic relation that connects if and 
Tf (Eqs.E]and[9]): 



if = - 2Z*_i + 2 7 k oS Tf + ^ ( CD,to - rf 



/■'oil 



ft'™ 



7 y 



CD. to 



Z.-D 



in 

'i-l 



exp 



on rpC 



V 



(D.4) 



Solving this equation for Tf gives a function Tf = 
g(lf^i,lf) that involves a LambertW function and can 
therefore not be expressed in terms of elementary func- 
tions. For correct normalisation, the new conditional 
probability distribution also has to be multiplied by the 



derivative 



The integration limits are the minimal and maximal 
values of l 7 f_ 1 , such that the order of events in Fig. [2] 
(i.e. regular oscillations) is still fulfilled for a given if. 
The upper limit turns out to be l mSLX , whereas the lower 
limit can be found as the solution to the polymer growth 
equation with Tf = Tfti- One finds 



jV 



fc ffLambertW 



1 



lVk oS 

(-7VcD,to*S. + 7™off + AOi) 
'c D: t fc° n - k oe 

toff 



• exp 



CD, to 



oir 



(D.5) 



Putting Eq. ID. II and Eq. ID. 31 together, one obtains a 
map from one probability distribution into the next one, 
half a period later: 



1 



7^off 



id jn 

T f = !i_JL\i< 



jk 



off 



dlf * cap 



P cap {T c = g{ltxMW-i) 



(D.6) 



The integrals can be computed numerically [38( and by 
iterating Eq. ID.6( a steady state probability distribution 
can be obtained. A typical starting distribution for the 
iteration would be a delta distribution or a uniform dis- 
tribution and for typical parameter values, a steady state 
distribution is reached after 3 to 5 iterations. 

This approach is valid, as long as the vast majority 
of events fall into the scheme depicted in Fig. i.e. as 
long as tf > and tf > t™ +1 (regular oscillations - see 
SubsecHSSJ. 

A more relevant quantity than l n , the length of a poly- 
mer at the time of nucleation on the other side, is the 
capping length l c , i.e. the amplitude of the oscillation. 
To obtain the probability distribution for this length, one 
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needs to perform another integration: 



P C (Z C ) = p(l c \l n )Pn(l n ) dl n . (D.7) 

The conditional probability distribution can be obtained 
from Eq. IC.ll when replacing T c with the appropriate 
function of l c and l n that can be derived from a similar 
relation as Eq. ID. 41 We obtain 



dT c (l c ,l n ) 



P cap (T c = T c (l c ,l n )\l n ). (D. 



The lower limit of the integration is the same as in the 
previous integration (Eq. ID.5[) . For the upper limit one 
needs to find l r ^_ l such that if = 0. This can be found by 



solving Eq. El with T c = 



this solution is larger than Z max , i max 
limit. 



- l c ) for l n . In case 
(Eq. 3]) is the upper 



APPENDIX E: PROBABILITY DISTRIBUTION 
FUNCTION FOR THE PERIOD IN THE 
OSCILLATORY STATE 

Using the iterative approach described in the previous 
appendix, the steady state probability distribution func- 
tion for the lengths of D-polymers at a given time point 
during regular oscillations can be found. Starting from 
this distribution, we will compute the probability distri- 
bution function for the period of the oscillations in this 
appendix. 

From Fig. [2] half a period T h (the time between two 
consecutive nucleations on opposing sides) is found to be 
T h = T n + T j and tnerefore p h (T h ) = J p{Tl = T h - 

T n \T n )p(T n ) dT n . T n can be replaced by T n = -±^l n 

and an additional probability distribution p(l d \l n ) can be 
introduced in order to obtain an equation involving only 
distributions that have already been calculated: 



are the smaller of { r yk «T h , Z max } for Z™ ax and Z^ ax is 

iv ( fe off ' 



the equivalent of Eq. ID. 51 I 



d 

max 



CD, to — JC 
( AC in\ 



For the probability distribution function of the period 
T = 2T h , one more integration is needed: 

P T (T) = [ P h (T h = T h ')P h {T h = T - T h ') dT h ' 
Jo 



(E.2) 



APPENDIX F: A LIMITING DETERMINISTIC 
MAP FOR THE CASE OF HIGH 
COOPERATIVITY 



For the case of deterministic nucleation and stochas- 
tic switching, we want to use the deterministic map as 
shown in Fig. [3] as a rough guideline for the dynamics 
of the system. In order to find this underlying map, we 
need to find a replacement for the parameter CE,th in the 
analytical expression of the map (Eqs. I A. If - IA.3|) since 
the MinE nucleation threshold does not have a meaning 
if we assume stochastic capping. A natural choice is the 
E-concentration at the median of the capping cumulative 
distribution function F cap (t). In the case of determinis- 
tic switching this is a step function in ce with the step 
at CE,th- For high rt cap , F cap (t) it is a steep sigmoidal 
curve. For the vast majority of cases (fc cap < — 



and fccap > 



ln(2) 



"E,0 



t d,E + 



™cap . d E 
L E,0 t 
-1 

), the 



median capping time tu lies within t d:E < tja < 
t d . With the integration of Eq. [CT] (and l n = 
■^X"(cD,to — c D,th)) one finds the median time to be = 



d,E 



*(- 



CE,0 



'E,0 



A(r 



cap J-d.Fj 



T, '"cap/, 

" / cap L 'E,0 



ln(2) 



which then leads to 



Ph{T h ) 



Pnuc(T f =T h 



1 



-l n \l d ) 



!i Jo ikott 

■ p{l d \l n ) dl d P n {l n ) dl n . (E.l) 



Pmic{T f \l d ) is given in Eq. E2 how to obtain p(l d \l n ) 
is described in the preceding appendix and a steady 
state distribution for P n (l n ) is also derived there (iter- 
ating Eq. ID.6p . The upper limits of the two integrals 



ce(£m) 



A{n cap + 1) 



'E,0 



^cap 



fccapCg Q V t l n (2) 



(F.l) 



We use this E-concentration as a replacement of CE.th 



in Eqs. lA.lllA.3l and plot the resulting maps in Fig. 10(d) 
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